An optical kaleidoscope using a single atom 
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A new method to track the motion of a single particle in the field of a high-finesse optical resonator 
is described. It exploits near-degenerate higher-order Gaussian cavity modes, whose symmetry is 
broken by the phase shift on the light induced by the particle. Observation of the spatial intensity 
distribution behind the cavity allows direct determination of the particle's position with approxi- 
mately wavelength accuracy. This is demonstrated by generating a realistic atomic trajectory using 
a semiclassical simulation including friction and diffusion and comparing it to the reconstructed 
path. The path reconstruction itself requires no knowledge about the forces on the particle. 



In a variety of pioneering experiments in the past few 
years |l] ||| it has been demonstrated and widely exploited 
that a single near-resonant atom can significantly influ- 
ence the field dynamics in a microscopic high-finesse op- 
tical resonator. Vice versa, the light field also influences 
the motion of a cold atom, which leads to an intricate 
dynamical interplay of atomic motion and field dynam- 
ics [f|^]. As a striking example, trapping of a single atom 
in the field of a single photon has become feasible Hf/J. 
This was experimentally substantiated by analyzing the 
characteristics of the measured output field. The time 
variation of the transmitted intensity shows very good 
agreement with theoretical simulations || of the con- 
fined three-dimensional motion of the atom in the cavity 
light field including friction and diffusion [9|-pd|. Carry- 
ing this analysis further it was even possible to associate 
piecewise reconstructed trajectories with recorded time- 
dependent intensity curves , utilizing the knowledge 
of the near-conservative potential. However, the recon- 
struction was possible only for atoms with sufficiently 
large and conserved angular momentum around the cav- 
ity axis and can only be done up to an overall angle and 
the direction of rotation. The reason was that only a 
single cavity mode, the TEM o mode, was used. Con- 
sequently, only a single spatial degree of freedom of the 
atom could be extracted directly from a measurement of 
the transmitted field. 

In this letter, we propose a new method to obtain two- 
dimensional information on the atom's position. We 
consider higher-order frequency-degenerate transversal 
modes. Typical examples are the Hermite-Gaussian 
(HG) or the Laguerre-Gaussian (LG) modes, which show 
in the transverse plane a rectangular matrix of intensity 
minima and maxima or a pattern of concentric rings, 
respectively. In such a configuration, the atom can redis- 
tribute photons from one mode to the other. Moreover, 
it induces mode frequency shifts and losses strongly de- 
pendent on the atomic position. In particular, the sym- 
metry of the intracavity field determined by the cavity 
and pump geometry will be perturbed and characteristic 
optical patterns containing information on the atomic po- 



sition appear. The effect reminds of a toy-kaleidoscope in 
which small objects in a symmetric arrangement of mir- 
rors create beautiful patterns. Our technique yields much 
more information on the atomic position and motion as 
compared to the single-mode case, and allows to extract 
the atomic position directly from a field measurement. 
To treat this problem quantitatively we generalize pre- 
vious semiclassical models of dynamical cavity QED in 
the strong-coupling regime to include finite sets of nearly 
degenerate eigenmodes. For a weakly saturated atom we 
derive a coupled set of equations for the mode amplitudes 
and the atomic center-of-mass motion. To be specific, 
let us consider a single two-level atom with transition 
frequency uj a and line width T (half width at half maxi- 
mum) moving inside a high-finesse optical resonator with 
transversal LG eigenmodes u pm (r), where p is the radial 
mode index and m is the azimuthal mode index lol , 



u pm (p, 9, z) = C pm cos(fcz)e w o : 
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where is the generalized Laguerre polynomial. The 
normalisation parameters, C pm , are chosen such that 
/ \u pm (p,0,z)\ 2 dV = dw 2 7r/4 = V 00 (the TEM 00 mode 
volume), where wo is the cavity waist and d the cav- 
ity length. At each spatial point the local atom-mode 
couplings are <? plTl (r) = go%> m (r), where go is the max- 
imum coupling of the TEMoo mode. For simplicity, we 
assume that the mirrors are ideally spherical and have 
a uniform coating. Then, all these modes have a com- 
mon eigenfrequency w and field decay rate k. However, 
the model can be extended to incorporate non-degenerate 
modes in a straightforward manner. The cavity is also 
assumed short compared to the Rayleigh length of the 
mode so that the wavefronts are approximately plane 
with z-dependence cos(27rz/A). The resonator is exter- 
nally driven by a coherent pump field of frequency lo p 
which pumps the modes with strengths r\ prn . Assuming 
low atomic saturation, we can adiabatically eliminate the 
internal atomic dynamics and treat the atom as a linearly 
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polarizable particle, which induces a spatially dep endent 
phase shift and loss. In the semiclassical limit VM, where 
we consider the center-of-mass motion of the atom classi- 
cally, we can derive the following set of coupled differen- 
tial equations for the mode amplitudes au (t) , the atomic 
position r a (t) and momentum p (t): 
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Here, M is the atomic mass, Uq — g^/{uj v — u> a ) is the 
single-photon optical light shift, 7 = Tg^/ (uj p — u> a ) 2 the 
spontaneous emission rate for a single-photon field, A = 
Lu p — lo the cavity-pump detuning, and x an d £pm are 
Gaussian random variables which model momentum and 
cavity field fluctuations, respectively. 
As an example, we consider an atom at rest at a given 
position r a coupled to the three degenerate cavity modes 
with (p, m) = (1, 0), (0, -2), (0, 2). It is then possible to 
solve equations @ for the mean stationary field ampli- 
tudes Q!p^(r a ), which in a parametric way depend on the 
atomic position. We get: 

s+at Vpm , ""I - 7 * / \r> / \ /o\ 

a p™ = iA^ + -K^ u p m{rMva) ' (3) 

where £q is the electric field at the position of the atom, 

0{ra) ~ (iA- K )-(i[/ + 7)E P , m Km(r a )| 2 - [ ' 

Note that the set of mode amplitudes u pm (r a ) at the po- 
sition of the atom may be written as a vector, which can 
then be rotated into the form (u e s(r a ), 0, 0). Since any 
linear combination of the modes can be considered as a 
mode as well, the atom-field dynamics for an atom at rest 
thus reduces to the case of a single mode u e g. This al- 
lows to derive simple analytical expressions for the steady 
state. For a moving atom, this approach must be gen- 
eralized, but still helps to find analytical expressions for 
friction and diffusion coefficients for the atomic motion 

Figure 1 shows the steady state field intensities for the 
empty cavity and for two different atomic positions. For 
the chosen parameters (rubidium atoms, (go, F, k) = 2t:x 
(16, 3, 1.5) MHz, (?ft0,%-2,?702) = 2n x (6.4, 0, 0) MHz, 
A = -2tt x 2.25 MHz and uj p - u a = -2tt x 114 MHz, 
leading to Uq = A, 7 — 2ir x 60 kHz) the atom dis- 
tributes photons between the cavity modes and changes 
their relative phases in such a way that a local maximum 
of the total field intensity is created near the position 



of the atom. By a change of the detuning between the 
pump laser, the cavity modes, and the atomic transition, 
a local minimum can also be achieved. Figure 1 shows 
the effect of the atom position on the shape and overall 
intensity of the stationary field distribution. 
This dependence of the cavity field on the atomic posi- 
tion suggests that measuring the cavity output field dis- 
tribution yields ample information on the atomic motion. 
In fact, it can be shown that the functions ap^(r a ) can 
be inverted almost everywhere to yield the atomic posi- 
tion in three dimensions. Of course one is limited by the 
common symmetries of all modes. For instance, for the 
system considered in Fig. 1, a 180° rotation around the 
cavity axis forms a symmetry operation. Hence, the re- 
construction of the atomic position from the cavity field 
will always yield two equivalent positions. Other sym- 
metry operations are a shift of A/2 in direction of the 
cavity axis and a reflection at the nodes or antinodes of 
the standing wave. Another limitation is that an atom 
cannot be detected close to the (transversal) nodes of 
the pumped mode (see Fig. la). However, it is possi- 
ble to determine directly from the photodetector signals 
whether the atomic position can be obtained or not: re- 
construction is possible if the transmission signal with an 
atom differs from the signal of an empty cavity. Also, the 
nodal areas are small and one is free to alternate rapidly 
between different pump geometries. Alternatively, one 
can change the pump geometry online when the atom 
approaches the nodal area of the pumped mode. 
Although in principle a full three-dimensional atomic 
trajectory can be reconstructed, the method encounters 
some complications. The longitudinal motion is in gen- 
eral too fast to be resolved experimentally. This amounts 
to replacing the coupling constant g$ by its longitudinally 
averaged value. A two-dimensional reconstruction still 
works in this case, even if the precise factor by which 
the coupling is reduced is unknown. Second, a single 
atom must redistribute enough photons among the cav- 
ity modes. This requires values of Uq of the order of the 
cavity decay rate k and hence the strong-coupling regime 
of cavity QED. Third, our arguments above are based on 
a stationary cavity field. For an atom moving in the xy 
plane, we thus have to assume that the cavity field fol- 
lows the atomic motion adiabatically. This implies slow 
atomic motion and not too large optical forces in this 
plane, which together with the constraint of strong cou- 
pling is tantamount to relatively low intracavity photon 
numbers. Fourth, in an actual experiment the exact in- 
tracavity photon number can only be deduced from the 
number of photons emitted by the cavity. The measured 
cavity output is subject to statistical fluctuations (shot 
noise), which become significant in the few photon limit. 
This limits the accuracy to which the cavity field, and 
therefore the atomic position, can be known. 
In the following we will demonstrate with a sample tra- 
jectory that all of these conditions can be met simul- 
taneously and that the numerical reconstruction of the 
atomic path is possible using existing high-finesse optical 



2 



resonators |6|j7|]. For simplicity we will assume a quasi 
two-dimensional situation where the atom is trapped in 
the longitudinal cavity direction at an antinode (z = 0) 
of the standing wave during the full interaction time. 
In a first step we create a sample trajectory for a single 
atom traversing the resonator by integrating the stochas- 
tic equations of motion (||) for given initial atomic posi- 
tion and velocity. This procedure includes all reactive 
and dissipative optical forces which the cavity field im- 
poses on the atom [Q, the back-action of the atom on 
the cavity field, as well as the momentum and cavity field 
diffusion. A resulting trajectory is depicted by the solid 
curve in Fig. 2. The atom enters the resonator from be- 
low. By chance, the atom encircles the cavity axis a few 
times before it is ejected again. 

The generated trajectory allows us to simulate a realistic 
cavity output signal. We assume an arrangement of 16 
photodetectors at the cavity output port each counting 
the numbers of photons detected in equally sized sectors 
extending from the axis to the edges of the cavity and 
covering an angle of 22.5 degrees. Due to the symmetry 
of the system, the signals from opposing detectors can be 
added without loss of information. To reduce the shot 
noise, we integrate the measured photon flux at each de- 
tector over a time interval of 100 cavity decay times 1/k. 
In Fig. 3 we plot the simulated photon flux for two out 
of the eight pairs of photodetectors. 
We can use the generated time-dependent field pattern to 
numerically reconstruct the particular atomic trajectory. 
This involves several steps. First, for each time step we 
determine the atomic position by a least-square compar- 
ison with a list of precalculated detector outputs corre- 
sponding to the steady-state field distribution a^(r a ) 
of any atomic position on a discrete spatial grid. Be- 
cause of the twofold spatial symmetry, we always obtain 
two equivalent points in that way. Hence, from a chosen 
initial point we have selected the one point out of these 
two that forms a continuous curve as a function of time. 
The spatial points obtained in that way are indicated by 
the crosses in Fig. 2. The initial position of the atom is 
usually known in experiments where laser-cooled atoms 
are injected into the cavity. For the pump geometry cho- 
sen here, there is a dark ring where the atom does not 
couple to the pumped mode. Close to this ring, the re- 
construction is difficult. As mentioned before, there are 
ways around this problem, but here the corresponding 
crosses are simply left out. 

According to the temporal noise of the measured photon 
fluxes (Fig. 3), the reconstructed atomic positions show 
a certain spatial spread. Since for the given parameters 
the stochastic forces in Eqs. (||) are much smaller than 
the Hamiltonian forces, we finally fit a smooth curve to 
this discrete set of data. The resulting reconstructed tra- 
jectory is shown by the dashed curve in Fig. 2. Note that 
the curve rotated by 180° forms a completely equivalent 
solution to the problem. Comparing the reconstructed 
with the original trajectory we note that for the depicted 
area close to the cavity axis the reconstruction works well 



with an accuracy which is of the order of an optical wave- 
length. 

The proposed detector arrangement might be con- 
structed by segment mirrors imaging onto an array of 
single-photon counting detectors. Probably more prac- 
tically is direct imaging on a high-sensitive high-speed 
camera. The construction of a suited cavity will be 
challenging. Scatter, misalignment and deformation of 
the high-reflectivity mirrors must be kept to a minimum 
to prevent breaking of the cylindrical symmetry, which 
could lift the frequency degeneracy of the modes by a too 
large amount. 

It is interesting to note that positive or negative values for 
to in the Laguerre-Gaussian modes correspond to actual 
left- and right-handed rotation of the light around the 
cavity axis. Hence, there is a striking similarity between 
Fabry-Perot cavities sustaining higher-order Laguerre- 
Gaussian modes and ring cavities. Similar dynamics as 
predicted in Ref. |l(| can therefore be expected in the 
system proposed here. 

An extension of the idea presented here is to use a cavity 
where modes with different longitudinal mode index are 
degenerate. In that case, one can choose a combination 
of modes with opposite parity to break the 180° rotation 
symmetry An example is a LG mode with even to de- 
generate with another one with odd to. In the extreme 
case of a confocal cavity, this will lead to a single field 
maximum near the position of the atom. 
In conclusion, we have shown how a high-finesse micro- 
cavity can be used as a real-time single-particle detec- 
tor with high spatial resolution. In contrast to conven- 
tional single-atom detection schemes, the cavity works 
as a phase-contrast microscope enhanced by the inherent 
multi-path interference of the high-finesse cavity. The 
method does not rely on fluorescence and, hence, can 
also work for particles without a closed optical transition. 
In contrast to the single-mode case, the information in 
the field pattern of sufficiently high-order modes should 
also allow one to spatially resolve several particles. It 
can be expected that the scheme presented here will first 
be demonstrated with single atoms in presently available 
high-finesse cavities, but applications extend beyond this 
system. The method should, in principle, also be appli- 
cable to large (bio) molecules in vacuum or even in pure 
watery solution. 
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FIG. 1. Transversal spatial intensity pattern (arbitrary 
units) of the stationary cavity field for the empty cavity (a) 
and two atomic positions indicated by the thick vertical line 
and a cross (b), (c). 
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FIG. 2. Central part of Fig. 1(a) with the simulated 
atomic trajectory (solid curve), reconstructed atomic posi- 
tions (crosses), and fitted atomic path (dashed curve). The 
atom enters with a velocity of 12cm/s and the total trajec- 
tory takes 5300 k -1 = 0.56 ms. The dashed circle indicates 
the dark ring of the pumped cavity mode uio. The waist of 
the TEMoo mode is wo = 29 fim. 
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FIG. 3. Simulated photon flux (number of photons iVph 
per 100 k _1 ) measured at two out of 8 photodetector pairs 
as a function of time (units of k _1 ). The inset depicts the 
arrangement of the detectors on the cavity output, shaded 
areas correspond to the plotted curves. 
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